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I. 


INTRODUCTION 


Modem  tactical  warfare  increasingly  uses  electronic  equipment  of  advanced  design. 
Although  radiation  threat  levels  at  the  tactical  level  are  relatively  low  compared  to  those  at 
which  strategic  or  space  electronic  systems  must  operate,  the  use  of  semiconductor  components 
requires  that  nuclear  vulnerability  be  assessed  and  the  equipment  hardened,  if  necessary. 
Statistical  methods  have  been  developed  to  provide  rapid  susceptibility  assessments  for  circuits 
containing  bipolar  transistors.  These  techniques  use  either  neutron-effects  or  manufacturer's 
data  to  predict  probability  of  survival  as  a function  of  neutron  fluence.  However,  these  methods, 
which  were  developed  primarily  for  assessment  of  new  equipment,  tend  to  include  large  errors  of 
the  second  kind;  that  is,  they  will  often  predict  that  a circuit  will  fail  where  in  reality  it  would  not 
fail.  In  order  to  use  these  methods  for  vulnerability  assessments  whose  goal  is  to  determine 
expected  values  for  circuit  failure,  it  is  necessary  to  extend  the  models  to  account  for  the 
distributions  of  the  initial  device  parameters.  Additionally,  the  current  models  are  applicable 
only  to  the  case  where  one  transistor  determines  circuit  performance.  When  circuit  performance 
is  determined  by  a critical  transistor  combination,  a Monte  Carlo  analysis,  generally  using  a 
network-analysis  code,  is  presently  required  to  determine  circuit  performance. 


Two  statistical  models  have  been  formulated'  to  predict  neutron  damage  to  a single  bipolar 
transistor:  one,  the  Electrical  Parameters  Model,  requires  only  the  electrical  parameters  given  in 
the  manufacturer's  specifications,  while  the  other,  the  Neutron-Effects  Data  Model,  makes  use 
of  the  radiation-effects  data  base  developed  by  the  Defense  Nuclear  Agency  (DNA)  and  Harry 
Diamond  Laboratories  (HDD.  Both  are  derived  from  the  Messenger-Spratt  equation." 


A,  ^ 

(1) 

1 K 

= + “T 

A,  U 

(2) 

where 

/3„  = initial  current  gain  (before  irradiation), 

= final  current  gain  (after  irradiation), 

(<)  = neutron  fluence, 

C = neutron  damage  factor, 

fr  = transistor  gain-bandwidth  product,  and 

K = fj-C  = normalized  neutron  damage  constant. 

The  Electrical  Parameters  Model  is  applicable  to  bipolar  transistors  in  general;  the  formula 
for  the  probability  of  transistor  failure  is  derived  from  equation  (2),  with  K treated  as  a lognormal 
random  variable  and  <t>,  fr  = frnwn  (minimum  specified  value  for  fr),  and  /3o  = jSmm  (minimum 
specified  current  gain  at  the  peak  of  the  gain  versus  current  curve)  as  constants. 

The  Neutron- Effects  Data  Model  is  appropriate  when  neutron-irradiation  test  data  are 


' D.  L.  Durgin.  D.  R.  Alexander,  and  R.  N.  Randall,  Hardening  Options  for  Neutron  Effects,  Harry  Diamond 
Laboratories  CR-74-052-1  tl5  November  1976). 

’ G.  C.  Messenger  and  J.  P.  Spratt.  The  Effects  of  Neutron  Irradiation  on  Germanium  and  Silicon,  Proc.  IRE,  4(t 
(1958),  1938. 
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(3) 


available.  It  is  based  on  the  rearranged  expression 


I _ >t/8o 
{A  + (t>)C  A + 4)' 


where  A - (Jj/fif,)/K  and  </)are  treated  as  constants  and  C as  a lognormal  random  variable.  Since 


A 


I 

CA, 


(4) 


or 


A.  . 

the  probability  of  failure  can  be  expressed  in  terms  of  the  distribution  of  the  random  variable  /So, 
which  is  also  lognormal  since  In/So  = -InC  — InA  is  normal.  Note  that,  in  contrast  to  the 
Electrical  Parameters  Model,  where  /So  and  K (or  C)  may  be  treated  as  independent,  the 
requirement  that  A be  constant  in  the  Neutron-Effects  Data  Model  leads  to  equation  (4),  which 
relates  A)  and  C.  This  equation  then  determines  the  relationship  between  the  probability 
distributions  of  Ai  and  C. 


The  Electrical  Parameters  Model  was  found  to  be  extremely  conservative  in  its  prediction  of 
failure,  mainly  because  it  treats  Ai  = /Smm  and  fr  = ^n„„  as  constants.  This  approach  is  helpful 
when  designing  or  hardening  equipment;  however,  in  a vulnerability  analysis  it  leads  to  large 
errors  of  the  second  kind  (i.e.,  the  circuit  is  predicted  to  be  soft  when  it  is  in  fact  hard).  Also, 
both  models  handle  two-transistor  combinations  through  a Monte  Carlo  approach  only. 

In  the  present  work  these  models  for  transistor-failure  probability  are  extended  to  include 
probability  distributions  for  the  initial  current  gains  and  to  allow  nonzero  origins  for  all  random 
variables  concerned.  In  fact,  the  formulas  developed  can  be  readily  adapted  to  any  type  of  input 
distributions.  Further,  these  models  are  generalized  to  consider  two-transistor  combinations. 
Test  cases  are  calculated  to  compare  the  failure  probability  curves  generated  by  these^  models 
with  previous  results. 


2.  ELECTRICAL  PARAMEIERS  MODEL,  WITH  p„  AS  RANDOM  VARIABLE 
Consider  the  expression 


_l_ 

A* 


(5) 


where  4>  is  constant,  and  and  C are  lognormal  random  variables  as  shown  in  Figures  1 and  2. 
(Unless  otherwise  specified,  all  constants  are  assumed  to  be  positive.)  Now  /So  and  C are 
assumed  to  be  independent.  If  this  is  not  the  case,  the  joint  density  function  of  Ao  and  C would 
be  needed:  with  this  the  ensuing  analysis  would  remain  valid.  Another  approach  is  available  when 
there  is  a high  degree  of  correlation  between  l/Ai  and  C. ' 


Standard  techniques  for  finding  the  distribution  of  a function  of  two  random  variables®  will 


' D.  L.  Diirgln.  D.  R.  Alexander,  and  R.  N.  Randall.  Hardening  Options  for  Neutron  Hffecis.  Harry  Diamond 
LatMiralories  CR-74-0.S2-I  ( I.S  November  I97fil.  p.  111-49. 

■’  Paul  G.  Hoel.  Introduction  to  Mathematical  Statistics.  4lh  Kdilion.  John  Wiley  & Sons.  Inc.  ( 19711.  p.  249  f. 


g(C) 


exp(-[ln(A)  - Anln)  - 


Figure  1.  Lognormal  distribution  for  ;3„. 


Probability  density  function; 
1 


je(C)  = — 

V2^er,.(C  - Q,i„) 
for  C > 


exp(-[ln(C  - C„|„)  - 


Figure  2.  Lognormal  distribution  for  C. 


be  used  in  determining  the  probability  density  Rearranging  equation  (5)  gives 

ft  = 

I + 4>poC 

from  which  it  is  clear  that  is  a decreasing  function  of  C when  /So  is  held  fixed. 
Then  the  joint  density  of  /3o  and  /3*  is  given  by 


7 


This  gives  the  probability  density  of 


P(/3J  = ] 


otherwise 


where  the  lower  limit  of  integration  arises  because  the  inequality 


Hj 


is  equivalent  to 


I - ■ 

When  it  has  been  determined  that  a given  circuit  will  fail  if  the  gain  of  a single  transistor  is  less 
than  some  threshold  value,  fir-  the  probability  of  circuit  failure  for  a given  neutron  fluence  d>can 
be  expressed  as 

Py(<b)  = Prob  < Pr) 


= f 


I 

where  ^r<  l/(Cn,|„<|»)  represents  the  transistor  gain  threshold  for  circuit  failure.  (If  At  2 l/(C„i„d)) 

then  Py{(i>)  =1.)  (-7^  m.n  /, 

However,  in  the  HONE  Electrical  Parameters  Model,  the  normalized  variable  K = frC  is 
used  in  place  of  C with  fr  constant  (equal  to  but  this  change  is  readily  incorporated  into 

our  derivation.  Since  C is  a lognormal  variable,  K is  also,  since  ln(A:  - K„t„)  = ln(C  - Qm)  + 

1^  IS  normal,  where  . In  fact,  it  is  clear  that  if  the  probability  density  function  for 

K IS  written  as 

then  we  must  have 


and  so 


Pk  = Pr  + lnJ5-, 
= O’f  , 


Jt  \Jt, 


where  g is  the  probability  density  function  for  C.  The  probability  of  failure  is  then  given  by 


8 


with  as  the  threshold  gain.  This  same  result  may  be  obtained  by  applying  the 

preceding  techniques  to  the  expression 


. 

+ K4> 

A> 

with  /3o  and  K independent  lognormal  random  variables  and  </>  and  fr  constant. 


3.  NEUTRON-EFFECTS  DATA  MODEL— CASE  OF  TWO  TRANSISTORS 
The  following  formula  is  used  for  a single  transistor' 

1 


li^ 


(A  + <t>K  ’ 

where  A = (fT/^)/K  and  </>  are  constants.  Since 


(9) 


^U/  C'/3o  ' 


we  may  replace  C by  \l{Afif,)  in  equation  (9)  to  obtain 

APo 


A + </) 


(10) 


Assume  that  Po  is  a lognormal  variable  as  before.  Clearly,  is  an  increasing  function  of  /3o 
and  so  the  probability  density  of  /3*  may  be  obtained  from  that  of  /3o  as  follows. 


PiPJ  = 


/(/So) 


dp„ 

/(/So) 


A ^ 


A A-  <i) 


(^4*) 


for  /3*  s 


AP„ 


A+  <j> 


(ID 


' D,  L.  Durgin,  D.  R.  Alexander,  and  R.  N.  Randall.  Hardening  Options  for  Neutron  Effects — Final  Report.  Harry 
Diamond  Laboratories  CR-74-052-1  (15  November  1976). 
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Note  that 


which,  by  including  a possibly  nonzero  origin  for  the  distribution  of  /3„.  extends  the  relationship 
between  distribution  functions  used  by  Durgin  et  al.'  to  develop  the  Neutron-Effects  Data 
Model. 


The  HONE  model  treats  one  transistor  at  a time.  If  there  are  two  transistors  in  a circuit  with 
failure  probabilities  Pf^(d>)and  independence  of  the  transistors  is  assumed  and  the  circuit 

failure  probability  is  calculated  from 

PM)  = I - [I  - Py^WIl  - ((/))]. 

This  report  generalizes  the  methodology  of  the  Neutron-Effects  Data  Model  and  the 
Electrical  Parameters  Model  to  treat  the  case  depicted  in  figure  3 where  a combination  of  two 
transistors,  T,  with  gain  /3,  and  T.^  with  gain  ^2  ■ is  used  in  a series  voltage  regulator  circuit  to 
amplify  an  input  current,  /j  (from  the  sensing  circuit)  to  some  final  current  If  (the  current 
necessary  to  maintain  a constant  voltage  drop  across  R,).  If  the  gains  /S,  and  jSj  decrease  because 
of  exposure  to  neutron  radiation,  then  the  sensing  circuit  increases  A to  compensate  for  the  gain 
loss  so  that  If  continues  to  produce  the  desired  voltage.  However,  practical  circuit  limitations 
place  a maximum  on  the  value  of  /(  that  can  be  input  to  T, . This  circuit  limitation  establishes  a 


' D.  L.  Durgin.  P.  R.  Alexander,  and  R.  N.  Randall.  Hardening  Options  for  Neutron  Effects — Final  Report,  Harry 
Diamond  Laboratories  CR- 74-052*1  (15  November  1976). 
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lower  limit.  /3r.  for  the  gain  product  below  which  the  circuit  will  malfunction.  Since  we  are 
interested  in  the  operation  of  the  circuit  after  irradiation,  we  denote  and  /Sz  by  and  , 
respectively,  and  compute  the  probability  of  circuit  failure  Prob(/3*,/3*,,  < /St-). 

For  the  Neutron- Effects  Data  Model,  the  final  current  gain  is  then  given  by 

_ '^iA)i  Afi 

A,  + <f>  A.,  + <{)  (/t,  + <t>)(A^  + <^) 

where  /So,  and  /Soj  are  lognormal  .-andom  variables  with  probability  densities 

./|(Aoi>  = r z :exp(-[ln(Ao,  “ A,m.n)  “ (13) 

- A,  min) 

for  I niin  * 

./nlAo-i)  = -7^= j r -exp(-[ln(  Ao-2  - A2m,n)  “ (14) 

v27rfr„^^(A„2  - A2min) 

for  Ai2  ~ A2min  ' 

Assume  that  Am  and  A02  are  independent.  The  expression 

rt  _ '^i‘^2A)I  Ai2 

* (A,  + (/))(A2  + </>) 

evinces  A...  as  an  increasing  function  of  Aic  when  Am  is  held  fixed.  Computation  of  the  joint  densi'v 
of  Am  and  A^  yields 

//(Aii.A.)  I 


{A,  + F 


./i(A0,)/2(A02) 


/ '^|A;AoI  \ U + *) 

y(A,  -I-  <j!))(A.2  + </))/  I fto  = ' 

(A,  + </>)(A;  + <t>)  / (A|  + (/))(A;  -f  (/>)A<, 

>4.^Am  \ 


MA,  * 


>4.^Am  " 

Integrating  over  Am  gives  the  probability  density  of  A*: 


A,  A2A0 


Ml  ■!■<). )M,  -i^  ft ife 

(A,  + </))(Aa  -t-  (/))  f 


P(/3l^)  = 


if  A*  ^ 


A,  AgAi  min  At  min 
(A,  + </))(Aj  -t-  (/») ' 


otherwise. 


1 


The  probability  of  circuit  failure  for  a specified  neutron  fluence  </>  and  gain  threshold  ft-,  is  then 
given  by 


PM)=- 


(A,  + d))(A,  + <i))  f f 

7^ J I 

’ - U, +wKA.  nil 

■f  /-i  > A|  A-yfiy  min  Ai  min 

{A^  + + (t>) 

A -c  /J  ^ mint^lmin 

" +d.)(A+d,)- 


(A,  + d))(A,  + 

a,a.A» 


4.  ELECTRICAL  PARAMETERS  MODEI^-CASE  OF  TWO  TRANSISTORS 
4.1  Constant  Initial  Gains 

Suppose  two  transistors  are  combined  as  in  the  previous  section  with  post-irradiation  current 
gains  given  by 

— = -^  + C,ih. 

A..,  M.. 

Mv,  A,.  ' 

where  C,  and  C,  are  independent  lognormal  random  variables  with  probability  densities 

I?, (C, ) = — ^ -exp(-[ln(C',  - - /[!,.  >72(r,?  ) for  C,  (17) 

\ 2n-<r,,,(C,  - 

gjtCj)  = — == — -exp(-[ln(C  - - /j.  ]V2<r,? ) for  Cj  > (18) 

\ 27r(r,/C2  - C,„„n) 

and  /3o,  = A. mm  and  A<«  = Aamm  are  constants.  Then 

A.^  = A^A.v, 

_ .^01  Ai3 

1 + Aih<^^’i  + An2<^>Q  + AoiAoz'i'^CQ 

is  clearly  a decreasing  function  of  C'j  when  C,  is  held  fixed.  The  joint  density  of  C,  and  A* 


computed  to  be 


MC,,/3J  = 


«i(C,)&(Q) 


3/3* 


j_r  . Mb ,] 


( /3oi  \ 

yl  + /3oid>Ci/  (1  + 


f =— !—  [ — il 


A)i  tr'  1 ' r A)i  A)2 .1^ 

The  probability  density  of  /3*  can  then  be  expressed  as 


P(/3*) 


^01  r^“ 

J c. 


I r /a.,/)..  .1 


I 


1 + 

A)i  A)2 


rg.(C,)ft  +A..«C,)  ')] 


If  0 < /3*  < ^ A,.<<>C, + A.2«^>Q™,„)  ’ 

0 otherwise. 


(20) 


The  probability  of  circuit  failure  in  this  case  then  becomes 


I r .] 

/3oi  I J 1 . ( 1 f ^01^02 ,1\  if,  in 

TJoMI  f-.  > + \Po2<t>  [pAi  +/3o,d>C,)  ' J j 


t'rW  = 1 


if  /8r^ 
1 if  ^ > 


PolPo', 


(I  ■*■  ^01  mInHl  ■*■  ^02‘^f'2mln) 
PoiPm 


(21) 


(I  + /3oid>C,  „,!„)(!  + ^02‘(>Qmln)  ’ 
for  a specified  neutron  fluence  (f>  and  threshold  gain  /Sj-. 

4.2  Variable  Initial  Gains 

Now  consider  the  case  where  is  as  expressed  in  equation  (19)  but,  in  addition  to  C,  and 
Cj,  j3oi  and  /3oj  are  lognormal  random  variables  as  well,  with  probability  densities  as  given  in 
equations  ( 13)  and  ( 14);  assume  that  all  four  variables  are  independent.  The  technique  presented 
by  Papoulis^  will  be  used  to  determine  the  probability  density  of  /3* . 


' A.  Papoulis,  Probability,  Random  Variables,  and  St(Kha.stic  Prwcsses,  McGraw-Hill  (l%5),  pp.  234-235. 
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1 


L. 


J 


The  transformation 


a _ A)|  A>2 

1 + + <<.a«Q  + 

C,  =C, 

Cj  = Cj 

A)i  ~ A)i 

has  associated  Jacobian 

flC,  flCj  <?/3oi  <5/302 

y=  1000 
0 10  0 
0 0 10 

^A)2 


_ ^01 

(1  + <^»3o,C,){I  + • 

Therefore,  the  joint  density  function  of  /3«,  C, . Q . and  floi  •*  given  by 


Then  pifi^)  is  computed  as  a marginal  distribution: 


P(^)  = i 


I 


Ci<t 

J c, 


8t(Q 


r !_1 

...J 

feCQ)  J p^c,<t> 


(1  4-  d<3o,C.) 


j3„,  1 - 


P,c,<i> 

1 - ^c.c,** 


1 - C,<t,p. 


ii:.  * ^'*). 


(22) 

dp^,dQdC, 


if  0 < /3*  s 


1 


min  ^2  min 

t 0 Otherwise. 

The  probability  of  circuit  failure  is  then  given  by 


r ' ' 

\Pt  mm  ^ 

Pt  mm 

gi(Q) 


J fee. 


1 -«.£■, 


■c.**  /3„,  1 C,<1)^ 


/i(A,.)/2 


I - C2<<>/3, 


“(i  * ^'*1 


(23) 

ld^„,dQdC,d/3^ 


for  a threshold  gain  ^ l/(C,minQmin*^*'‘)  (otherwise  - •)• 

5.  COMPUTATION  OF  MODEL  PARAMETERS 

In  general,  suppose  a sample  mean  iand  sample  standard  deviation  Sj.  have  been  computed 
for  a random  variable  x known  to  have  a standard  lognormal  distribution 


/(Jt)  = 


1 


V^o-Jf 


exp[-(lnJf  - for  x > 0. 


For  this  distribution. 


m s mean  = exp^/x  + ^ , 
s * standard  deviation  = + - |) 


= m Vf"'  - 1 
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To  estimate  ix  and  a,  set 


x=  m = exp^M  + yj . 
s,  = J = x Vf  - I . 

Then  we  get  the  following  formulas  from  which  to  compute  <r  and  fx: 

■ 

M = In  v - ^ 

= In  i - In  yjl  + ^ 

= lJ-y£=J\. 

WP  + slJ 


(24) 


(25) 


If  the  origin  of  the  lognormal  distribution  is  shifted  to  some  value,  say  x„,„ , then  x - x„,|„  has 
a standard  lognormal  distribution,  so  these  same  formulas  can  be  used  as  long  as  x is  replaced  by 
- ■’^n.in ; that  is. 


(7 


• Viii[ 


I + -r- 


(X  - X^,„Y 

(■f  - 


(26) 

(27) 


Equations  (26)  and  (27),  with  x — ^ and  x = C (or  K),  are  used  to  compute  the  parameters  of  the 
lognormal  distributions  of  /S#  and  C (or  K)  in  the  extension  of  the  Electrical  Parameters  Model  in 
section  2,  and  to  compute  Ma,,  < Mr, < and  in  the  two-transistor 

cases.  ‘ 


The  computation  of  A in  the  Neutron-Effects  Data  Model  requires  the  statistics  of  C as  well 
as  /3o.  Since  C = l/(/4)3„),  we  have  InC  = -ln/8o  - In4  and  £tlnC]  = -£[ln/3o]  - ln/1  as  the 
relation  between  the  corresponding  expected  values.  This  gives  A = = exp[— (£lln/3o]  + 

£tlnC])].  Thus,  A can  be  estimated  from 

A = expl-(ui^  + b^c)l  (28) 

This  formula  requires  estimation  from  the  data  of  the  means  of  ln/3o  and  InC  rather  than  of /So  and 
C.  However,  if  /So  and  C have  unshifted  lognormal  distributions 

^ — ;rexp[-(ln/3i)  - A«B.)V2(r2j  for  /a„  > 0 

V 

V(C)  = —= — -exp(-(lnC  - fXfYUaf]  for  C>  0 

V 2n(r,  C 
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A can  be  estimated  from  the  statistics  of  /So  and  C themselves.  For  it  is  clear  that 

t'[ln/3ol  = 

£[lnC]  = He 


so  that 


Therefore,  we  have 


Hr  = - InA 


A = exp(  -(flln/Soi  + fllnC))] 

= exp[-/xa„  - H,  ] 

exp(o-gJ 

exp(/itf„  + ia-l)eKp{H(  + i®’?) 


where  s is  the  standard  deviation  of  the  lognormal  distribution  of  /So  - Consequently,  in  this  case 
an  estimate  of  A may  be  made  from 


' [teV.,], 


where  /So  and  \ are  the  sample  mean  and  sample  standard  deviation  of  /So,  and  C is  the  sample 
mean  of  C.  This  formula  agrees,  except  in  the  notation  used,  with  that  given  by  Durgin  et  al.' 


6.  RESULTS  OF  CASE  STUDIES 

Test  cases  were  examined  for  each  of  the  model  extensions  derived.  The  numerical  values 
of  the  statistical  parameters  (tables  I and  II)  were  obtained  from  manufacturer’s  specifications  or 
from  the  DNA/HDL  radiation-effects  data  base.  The  following  cases  were  considered. 

1.  Extended  Electrical  Parameters  Model  applied  to  a single  transistor  (2NI485) — /So  and  C 
lognormal  random  variables. 

2.  Shifted  HONE  Electrical  Parameters  Model  derived  for  a two-transistor  combination 
(2N1183  and  2N697) — C,  and  Cj  lognormal  random  variables,  /Soi  = /3,niin  and  /S02  = /Samm- 

3.  Same  as  case  2,  but  with  a different  value  of  /3,niin-  This  shows  the  sensitivity  of  the 
HONE  Model  to  the  selection  of  /Smm  • 

4.  Neutron-Effects  Data  Model  derived  for  a two-transistor  combination  (2N1183  and 
2N697) — /3oi  and  /Soj  lognormal  variables.  A,  and  Aj  are  computed  from  C, , Cj , /S#, , and  /3o2  data. 

5.  Extended  Electrical  Parameters  Model  applied  to  a two-transistor  combination  (2N1183 
and  2N697) — /So, , /Soj , C, , and  Cj  all  lognormal  random  variables. 

' D.  L.  Durgin.  D.  R.  Alexander,  and  R.  N.  Randall.  Hardening  Options  for  Neutron  Effects — Final  Report,  Harry 
Diamond  Laboratories  CR-74-052-1  (15  November  1976). 
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TABLE  I 


Input  Parameters  for  Single-Transistor  Cases 


Case  1 

Extended  Electrical  Parameters  Model 
with  Pa  and  C variable 

HONE  Electrical 
Parameters  Model 

Shifted  HONE  Electrical 
Parameters  Model 

J3„  =49.1 

Po  ~ ^mln  ~ 35 

^0  ~ ^mln  “ 35 

6.10 

C = 2.42  X 10-’" 

C = 2.42  X 10"’" 

Anin  = 35 

s,.  = 3.77  X 10-'” 

s,.  = 3.77  X I0-'» 

C = 2.42  X 10-'^ 

= 0 

C„,„  = 1.289  X |0->" 

s,.  = 3.77  X |0-“ 

0 

II 

0 

a 

Qin  = C - 3v  = 1.289  X |0-‘^ 

0 

II 

Note;  Single  transistor — 2N1485;  <i>  is  in  n/cm*,  C in  cmVn,  and  is  dimensionless. 


TABLE  II 

Input  Parameters  for  Two-Transistor  Cases 

Case  2 Case  3 

Shifted  HONE  Electrical  Parameters  Model  Shifted  HONE  Electrical  Parameters  Model 
Two-Transistor  Combination  Two-Transistor  Combination 

C, . Cj  variable;  /Sm  , ^02  constant  C, , Cj  variable;  fto, , /3o2  constant 


T,  : 2N1183 

7i:  2N697 

Same  input  parameters  as  Case  2 

1 

$ 

except  Pi  n,i„  = 30  i 

C,  = 3.0  X 10-'^ 

Q = 1.8  X 10-’* 

1 

i 

1 

s, , = 4.0  X 10-'* 

= 2.3  X 10-'« 

min  ^1  “ 

^2  min  ~ ^2  ~ ^^(‘2 

1 

1 

1 

li 

= 1.8  X 10-'* 

= l.II  X 10  '* 

A»i  = Pi  min  ~ 20 

A)2  — Aamln  — 'W) 

666 

J 

1 

Case  4 

Case  5 

1 

» 

Neutron-Effects  Data  Model 

Extended  Electrical  Parameters  Model  ' 

\ 

Two-Transistor  Combination 

Two-Transistor  Combination 

All  < Ak  . <7, , Cj  variable  ’ 

T,;  2N1183 

Tj:  2N697 

T,;2N1183  T,;  2N697  j 

\nr^  = -33.495 

InT^  = -33.959 

Aoi  = 55.4  Pm  = 77.0  ‘ 

Aoi  = 55.4 

P02  — 77.0 

= >7.9  = 10.8  1 

= 17.9 

5^  = 10.8 

Almln  = 20  Almln  =40  \ 

Al  min  — 20 

ftmln  — ^ 

C,  =3.0X10-'*  Cj  = 1.8  X 10-'*  1 

= 3.%7 

ln/3oj  = 4.334 

V,  = 4.0  X 10-'*  V,  = 2.3  X 10-'*  1 

A,  = 6.666  X 10'* 

A,  = 7.3447  X 10'* 

C,„m  = 1.8  X 10-'*  Q„,n  = 1.11  X 10-'*  j 

666 

/3t  = 666 

Note;  </>  is  in  n/cm',  C,  and  Cj  in  cmVn,  /3oi  and  Pm  dimensionless,  and  A,  and  A,  in  n/cm*. 
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Figure  4 indicates  the  significance  of  the  model  extension  for  the  one-transistor  case  that  has 
been  developed  in  this  paper.  There  are  two  aspects  of  the  improved  one-transistor  model 
developed  here  which  contribute  to  this  difference:  (I)  /3o  is  a random  variable  rather  than  a 
constant,  and  (2)  the  lognormal  distribution  for  C is  allowed  to  have  an  origin  other  than  at  zero. 
Figure  4 indicates  that  the  major  contributor  to  this  difference  is  the  incorporation  of  the 
probability  distribution  of  because  the  middle  curve,  which  represents  the  results  obtained  by 
j the  shifted  HONE  Model  ()3„  constant),  is  practically  indistinguishable  from  the  curve  for  the 

unshifted  HONE  Model. 

The  probability-of-survival  curves  for  the  two-transistor  cases  are  shown  in  Figure  5.  The 
!.  curves  for  cases  2,  3,  and  4 were  obtained  by  numerically  integrating  the  appropriate  expression 

; for  Pf  (4>)  for  a sequence  of  values  of  <|).  However,  in  case  5,  Py((t>)  is  given  by  a four-fold 

integral — equation  (23) — which  does  not  lend  itself  to  computer  evaluation.  This  curve  was 
i:  developed  by  a Monte  Carlo  technique:  for  each  value  of  4>,  the  distributions  of  /3oi , /3o2 , C, , and 

, €-2  were  sampled  times;  for  each  sample,  was  calculated  from  equation  (19)  and  compared 

f to  j3r;  the  number  of  times  that  jS*  was  less  than  Pt  was  counted.  This  count  divided  by  N*  then 


NEUTRON  FLUENCE  # In/cm*) 


Figure  4.  Comparison  of  extended  model  with  HONE  Model — one  transistor. 
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NEUTRON  FLUENCE  ^(n/cm*) 

Figure  5.  Comparison  of  models  for  two-transistor  combination. 


provided  an  estimate  for  = Prob(/3*  < j9r).  This  was  done  for  a sequence  of  values  of  </>, 

with  N*  chosen  large  enough  (usually  A/*  = 2000  was  sufficient)  to  ensure  convergence  of  the 
estimate. 

The  position  of  the  curve  for  case  2 compared  to  that  for  case  3 indicates  the  importance  of 
proper  selection  of  the  input  parameter  )3imin  in  the  HONE  Model.  It  was  to  be  expected  that  the 
curves  for  cases  4 and  5 would  compare  favorably,  since  each  of  these  cases  includes  the 
statistics  of  all  four  random  variables  /3oi . fim . Q . and  Cj . 

7,  CONCLUSIONS 

The  conservatism  of  the  HONE  Electrical  Parameters  Model  for  the  probability  of  failure  of 
a single  transistor  due  to  neutron  fluence  has  been  reduced  by  including  the  probability 
distribution  of  the  initial  current  gain  /3o  and  allowing  for  a nonzero  origin  for  the  distribution  of 
the  neutron-damage  factor  C.  Test  cases  indicate  that  these  improvements  lead  to  a significant 
change  in  the  survival  probability  curve,  with  most  of  this  difference  due  to  the  inclusion  of  the 
probability  distribution  for  /So. 

The  HONE  model  for  probability  of  failure  due  to  neutron  fluence  has  been  extended  to 


j 

i 

I 
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consider  a circuit  whose  function  depends  on  the  product  of  two  transistor  gains.  The  following 
models  have  been  derived  for  this  two-transistor  combination. 

1.  Electrical  Parameters  Model,  with  neutron-damage  factors  variable  and  constant  initial 
current  gains 

2.  Extended  Electrical  Parameters  Model,  with  neutron-damage  factors  and  initial  current 
gains  all  variable 

3.  Neutron-Effects  Data  Model 

In  all  these  models  nonzero  origins  for  the  distributions  of  the  random  variables  concerned  are 
allowed. 

Finally,  formulas  have  been  derived  to  estimate  the  parameters  of  the  probability  distribu- 
tions from  transistor  test  data. 
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